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Abstract 

In this paper we present a mathematical model for estimate the col- 
lapse time of a gas bubble in a vane of a oil gerotor pump. This amount 
of time cannot be greater of the total time spent by the pump for filling 
and then emptying out a vane in a single revolution, otherwise there is a 
loss of lubrication between internal and external gears. We assume that 
oil is incompressible and viscous, the bubble has a spherical shape and it 
is not translating into the external fluid. The analytical treatment of the 
model shows that the Navier-Stokes equations for the velocity field of the 
oil can be reduced to a single non linear ordinary differential equation for 
the variation in time of the bubble radius. The collapse time estimated by 
a numerical resolution of this equation and the collapse time calculated 
from an analytical resolution of the linearized equation are substantially 
equal. 

keywords: cavitation, Navier-Stokes equations, spherical coordinates, 
numerical and analytical resolution. 

1 The problem 

We consider an oil gerotor-type pump for industrial applications, with an 
internal maximum pressure of about 15 bar. Experimental results show that 
for a rotational velocity of about 2000 rpm some phenomena of cavitation can 
occur. In this work we don't consider a description of the arise of cavitation 
(see e.g. PP), but we want to estimate the time, which we call collapse time 
tc, spent by a bubble of gas to reduce to zero its radius. In this way, we can 
deduce a maximum rotational velocity to avoid loss of lubrication due to the 
presence of bubbles in the geometrical zone where profiles of pignon and crown 
are coincident (see Figure P). 
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Figure 1: Pignon and crown profiles in gerotor pump. 



Assume that the gas (air or oil vapours) bubble has a spherical shape, with 
radius function of time: R — R[t). Therefore the mathematical problem is to 
compute the time tc such that R{tc) — 0. The bubble is compressed by the 
external oil, which in general has a pressure p greater than pressure Pg of the 
internal gas. For simplicity, we assume that near the bubble the velocity field 
of the oil is radial, that is 



where is the radial versor exiting from the centre of the bubble, r is the radial 
distance of the oil material point from the centre, and t is time variable (see 
Figure 2), and that the bubble is not translating into the fluid. Also, we don't 
consider external mass forces acting on the oil. 

Our aim is to use Navier-Stokes equations for the oil flow near the bubble, 
so that we can use the velocity v(R,t) at the bubble-boundary for computing 
the unknown function R{t). Usual expression of these equations in a three 
dimensional cartesian system and without external forces is (see e.g. 



where p is the density, /i the viscosity and p the pressure of the fluid. For 
symmetry reasons due to the geometrical conflguration, we'll see that this form 
is not the best for the treatment of the problem. 

2 Pressure of fluid near the bubble 

Before using Navier-Stokes equations, we try to found an expression for the 



V = v{r, t)er 



(1) 
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Figure 2: Radial velocity field near the bubble. 



pressure p of the fluid external to the bubble. For symmetry, we assume that 
near the bubble p is a function of the radial distance r. Let Rq the initial radius 
of the collapsing bubble. 

Let Pm the mean pressure, which we can assume to be constant, of the fluid in 
the pump. Then, for r = i? = 0, that is when the bubble is collapsed, p is equal 
to Pm, while for r = R = Rq p can be considered as the difference between pm 
and the pressure Pgo due to the gas contained on a sphere of radius Rq. Prom 
the perfect gas law we have 

3noRT 

where uq is the mimber of gas moles, R is the gas universal constant and T the 
gas temperature. If we assume that p is linear on r, the simplest expression for 
p to satisfy the previous conditions is 

R() — R , , 

p = cRr + — Pm (4) 

where c is a constant which can be computed imposing the condition;) = Pm—PgO 
when r = R = Rq. The result is 

(Prn_ 3noRr\ RO-R ..s 

Note that, as expected, p is a radial symmetry function. 
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3 Use of spherical coordinates system 



The physical nature of the problem and its mathematical modelization induce 
us to substitute the cartesian coordinates with a spherical coordinates system 
with origin at the centre of the bubble. These coordinates are the radial distance 
r from the origin, the azimuthal angle (longitude) a, with < a < 27r, and the 
elevation from the plane of the equator (latitude) 6*, with < 6* < tt (see |3]). 
We have to translate the Navier-Stokes equations Q from cartesian to new 
spherical coordinates in the case of radial velocity field v(x, t) = w(r, t)er and 
radial pressure p = p{r, t). Note that r'^ ^ x'^ + y'^ + z"^ , therefore 



d^r = -, dyr = -, d^r = - (6) 
r r r 



Also we have 



e. = - (7) 

r 

Then, denoted by Vi the i-ih cartesian component of the velocity v, from the 
chain rule the following formula for the k-ih. component of the gradient of Vi 
holds: 



dkV; = ^ 1 

or 



From these identities, the cartesian divergence V • v has the following shape in 
the new coordinates: 
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V-v^ Va.^;.^^ + -^ (9) 
^-^ or r 

i—i 

Now we apply the hypothesis of incompressibility for the fluid; from V • v = 
follows 

dv 2 

jr- = —V 10 

or r 

which can be resolved by variables separation, so that 

v{r,t) = l^ (11) 

with / arbitrary function of t. If we assume that the bubble boundary has 
velocity R{t) equals to the external fluid velocity field when r = R{t) (physically, 
this corresponds to the fact that fluid doesn't enter into the bubble deforming 
its geometrical shape), we have 

R{t) = viRit),t)^I^ (12) 

and therefore f{t) = R^{t)R{t). Then, from the expression for the velocity 
fleld is 
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, , R^{t)R{t) 
^{r,t) = ^^^^' er (13) 

Now we can write the partial derivative respect time variable, necessary in 
Navier-Stokes equations: 

^ = ^{R'i^ + ^R^')^r (14) 

Then we translate in spherical coordinates the cartesian expression v • Vv = 
(Sfe=i "^kdkVi) ■ From identities lO and from Vk = f(r, i)(e,.)fe = v^, we 
obtain 



dv X 



VkdkVi = -T^v^ + V 
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k , „,2 



(15) 



Summing on k and using x"^ + + = r"^ , e.^ = ^ , the following expression 
holds: 

„ dv ^x y z\ dv 
V • Vv = 



—y (T. 1 1) ^ _„e (16) 

dr xr'r'rJ dr 



dr \r'r'rJ dr 
From and (|10|l we can deduce the new form of the convective term in fluid 
dynamics equations: 



V Vv= ( -^i?*(t)i?2(t) ) e, (17) 



Now, using standard computations (or see (0]) for a calculation in the case of 
spherical symmetry), from © and pU|l we can deduce the expression for the 
laplacian of v in spherical coordinates: 

Av = 5] 9Lv -iQ^ + -Q^)er^ -,ver ^ -^R\t)R{t)er (18) 
fc=i ^ ^ 

Finally, we derive the form for the gradient of the pressure p: 

(Vp). = |^^=.fc-^')i?^ (19) 
^ ^ ' dr r \Rl AttRI J r ^ ^ 



so that 



„ ,'pm 3noRr\ 
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4 An ordinary differential equation 



In this section we derive an ordinary differential equation which translate 
Navier-Stokes system into a new form using radial symmetry and the spherical 
coordinates. 

We substitute the differential expressions in original fluid dynamics equations 
with their new forms computed in the previous section. For simplicity, let 
a = ^ - ^^1?, with a constant. Therefore, using lfn|) . l(T7|) . ifTFI) and in 
equations we can write 

P (^^!^±^ . ^E^E^^ = -aR + A^R^R (21) 

For r — R this equation describe the evolution of the bubble radius; simplifying 
similar terms and multiplying by i?^, we obtain the following equation for the 
unknown function R{t) 

pRR^ - 2^iR + aR^ ^0 (22) 

Note that the leading differential term RR^ is the same than that of Rayleigh- 
Plesset equation of bubble dynamics in cavitation (see ^). The previous non 
linear second order equation seems not analytically integrable with classical 
methods; probably an hard calculation using Lie symmetries technique might 
be useful for informations about expression of an exact solution (see e.g. (2|). 
But in the next sections we deduce results of industrial interest using numerical 
and approximate analytical resolution. 

We want to solve an initial value problem for H22|l with conditions 

i?(0) = Ro, R{0) = (23) 

The condition on derivative seems quite realistic. Note that a not null constant 
function R{t) = k Wt is not a solution of the equation, as expected for physical 
and technical reasons. 



5 Physical considerations 

For a computation of the coefficient a in equation H22|) the main step is the 
estimate of the number of moles. As first, but quite realistic, approximation we 
consider that the gas contained in the bubble is atmospheric air. Let Ro = 0.05 
cm the initial radius of the bubble, a value compatible with the dimensions of a 
body pump of industrial interest. Using the value pair — 0.01177 g/cm^ as air 
density and the value Wair — 28.97 g/mol as air molecular weight, we calculate 
the air mass Mair contained in a sphere of radius Rq, so we can know the number 
of moles no = j^""' (see e.g. |S] for details). The resulting value of a is about 
10^ dyne/cm'^. As fluid we use light oil with p — 8.2 g/cm^ and p = 0.0287 
St (at standard temperature of 300°/'i'). As mean pressure of fluid into vanes 
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we consider the value pm = lO'' dyne/cw? « 10 bar, that is the value usually 
reached at the halfway during the angular rotation of a single vane pump. 

6 Numerical resolution 

Consider now the numerical resolution of our differential problem. We have 
used the software Mathematical ver. 5.2, and its built-in function NDSolve, 
which applies adaptive finite differences schemes like Adams and Backwards 
Differentiation (see [7] for informations about these methods) with high degree 
of accuracy (see |S] for details). The result is shown in Figure 3. 
The collapse time tc is about 0.00141 sec. At a frequence of 2000 rpm, during 
an interval of time equal to tc the pignon rotates of about 17°, therefore if 
bubble is produced near the geometrical zone where pignon and crown have 
conciding profiles (right zone of Figure^, the probability of loss of lubrication 
and physical contact between rotating parts can be not null. 




0.0007 0.00141 



Figure 3; Profile for numerical solution; variable in horizontal axis is time in sec. 



7 Analytical resolution 

In this section we approach an analytical approximated resolution for equa- 
tion (|22|l . From previous numerical resolution the key information is that 
tc « 1, therefore an approximation based on some type of expansion respect 
the time variable and centered at t = is well justified. For this aim we consider 
an application of Implicit Function Theorem to the case of symbolic resolution 
of differential equations (see for details). 

From hypothesis we have i?(0) = i?o and i?(0) = 0, so from (|22|) follows that 

m = (24) 

p 

Consider the polynomial function F{t, R, R, R) — pRR^ — 2p,R + aR^ , where R, 
R and R are to be considered independent variables and not functions. From 
previous formula, we have 
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F(0,i?o,0,-^) = (25) 
Then we compute all the partial derivatives of F at point (0, i?o, 0, — ^^): 

dtF = 0, OrF ^ aRl, djjF ^ -2p, dj^F = pRl (26) 

As dj^F 7^ 0, from implicit function Theorem (see e.g. ^01) we can write the 
following approximated but explicit formula for the variable R 



R = ^ - -{R~ Rq) + ^{R-Q) (27) 



«-^-^(i?-i?o) + %, 

p p pRo 

from which we obtain a linear constant coefficients ordinary differential equation 
for the function R{t): 



pRlR - 2pR + aRlR = (28) 

The analytical solution of the initial value problem can be found with the clas- 
sical method of the characteristic equation, and its shape is quite complicated: 











e "o" Ro 


y/p^ - aR^p 


1 + e ''op 


- e ^0" ^ + /i 



2^p^ - aR^p 

(29) 

But numerical experiments testify that its second order Taylor expansion cen- 
tered aXt — Q gives a very accurate approximation, so that we can establish that 
the shape of the approximated analytical solution of the differential problem of 
bubble collapse is 



R{t) = i?(0) + R[Q){t - 0) + ^i?(0)(t - Qf = 

2 2p 



Rf) 



(30) 



The collapse time tc in this case is given by resolution of the algebraic equation 
R{tc) = 0: 



tr = \ — 



(31) 



Substituting the values of p and a, we found tc = 0.00128 sec, which matches 
the previous value estimated by numerical resolution until millisecond. Note 
that viscosity p, in this approximation, is not relevant for computing collapse 
time. If a Taylor expansion of third order is used for the analytical solution of 
(|2HI) , the coefficient of is 



as 



ap 

WRq 



(32) 
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and the numerical contribute of the cubic term to reduce R(t) is very small, 
last^l < 10-'^ cm. 

In Figure 4 we compare the two profiles of R{t) obtained by numerical and 
analytical methods. 




Figure 4: Light line is the analytical profile. 



8 Conclusions 

The analytical resolution gives a formula for estimate the collapse time, and 
in particular from (|31|l and H32I) , we can deduce that a fluid with a lower density, 
or with a greater viscosity, or a pump with a greater mean pressure in vanes can 
give benefits for a reduction of time permanence of bubbles. The approximated 
calculations in analytical case give results matching very well those of numeri- 
cal resolution. But in the second case the previous technical informations about 
density and pressure can be obtained only by an amount of computational sim- 
ulations. 

The mathematical model can be improved, in particular assuming for the bub- 
ble a geometrical shape more complicated than the spherical one. In this case 
the function /(t) of equations 1)11(1 and l|12(l should have a more complicated 
dependence on R{t), therefore the differential equation could become more dif- 
ficult to resolve. 

The non linearity of the differential equation H22|l induces to consider more sub- 
tle techniques for searching a non approximated analytical resolution. Probably 
the method of Lie symmetries generating a two dimensional Lie algebra could 
transform the original equation to an integrable form. 
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